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' Abstract Observations of the inner radian of the Galac- 
' : tic disk at very high energy (VHE) gamma-rays have re- 
I ' vealed at least 16 new sources. Besides shell type super- 
nova remnants, pulsar wind nebulae (PWN) appear to be 
a dominant source population in the catalogue of VHE 
gamma-ray sources. Except for the Crab nebula, the 
newly discovered PWN are resolved at VHE gamma-rays 
to be spatially extended (5-20 pc). Currently, at least 3 
middle aged (i > 10 kyrs) PWN (Vela X, G18.0-0.7, and 
G313.3-f0.6 in the "Kookaburra" region) and 1 young 
PWN MSH 15-5^ {t = 1.55 kyrs) have been identified 
to be VHE emitting PWN (sometimes called "TeV Ple- 
rions"). Two more candidate "TeV Plerions" have been 
identifed and have been reported at this conference 
In this contribution, the gamma-ray emission from Vela 
X is explained by a nucleonic component in the pulsar 
wind. The measured broad band spectral energy distri- 
bution is compared with the expected X-ray emission 
from primary and secondary electrons. The observed X- 
ray emission and TeV emission from the three middle 
aged PWN are compared with each other. 
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1 Introduction 

The discovery of extended gamma-ray emission from pul- 
sar wind nebulae (PWN) opens an exciting possibility 
to study the acceleration of particles in ultra-relativistic 
shocks. Acceleration at relativistic shocks is of relevance 
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for the understanding of the mostly non-thermal emis- 
sion seen from relativistic jets in active galactic nuclei 
(AGN) and gamma-ray bursts (GRBs). While the Lorentz 
factor assigned to the relativistic flow in AGN and GRBs 
is currently believed to be in the range of ~ 10 . . . 100 
based upon arguments of opacity against pair produc- 
tion, the Lorentz factor of pulsar winds 7 is largely un- 
determined. Using the potential drop of the last open 
field lines as a constraint: 7 < e<l>open/ {rrieC?) = 2 x 
10^(P_i3/P)^/^ with P indicating the rotational period 
and dP/dt = P_i3 x 10^ ^■^ s/s the rate of slowing down. 
The Lorentz factor has been determined in model depen- 
dent ways to be in the range of 10^ — 10^ for the Crab 
nebula (see e.g. |21 and [3]). Even if this value is pos- 
sibly smaller in other PWN, we are witnessing particle 
acceleration in ultrarelativistic shocks driven by the en- 
ergetic pulsar wind. The total energy available by slow- 
ing down of the pulsar (assuming a moment of inertia 
/ = /45IO45 g cm2): E = 10^4 h^P.^^/P^ erg/s. 
Theoretical calculations of Fermi- type acceleration in rel- 
ativistic shocks using widely different approaches (see 
e.g. 0], [S]) lead to the consistent conclusion, that (largely 
independent of the detailed parameters of the shock) par- 
ticles in the downstream plasma are accelerated to fol- 
low a universal power-law for large Lorentz factors of 
s — 2.2 . . . 2.3 which is compatible with e.g. the observed 
X-ray synchrotron spectrum of e.g. the Crab nebula. 
However, recent progress in the understanding of the 
downstream turbulence spectrum, a revision of this "uni- 
versal acceleration" picture may be necessary (see e.g. 

[Z|)- In the more recent calculations, the particle 
spectrum is found to be softer than in previous calcu- 
lations (see above), to deviate from a universal power- 
law, and to show a stronger dependence on the Lorentz 
factor of the upstream medium and on the particular im- 
plementation of downstream turbulence. Similar conclu- 
sions have been found independently with particle-in-ccU 
(PIC) simulations (0, 0, [El)- 

The presence of ions in the pulsar wind obviously compli- 
cate the structure of the shock and have been described 
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e.g. in [ni, lEl, |H|- Hoshino et al. (1992) [H] found that 
an admixture of ions in the wind can lead to accelera- 
tion of positrons in the downstream region by resonant 
absorption of magnetosonic waves emitted by the gyrat- 
ing ions. In a recent PIC simulation [Tn| a larger ratio 
of ion mass m.i to electron mass me in the simulation 
has been used (previously, PIC simulations had been 
limited to values of TOi/mg « 20) and various energy 
fraction of ions have been considered. The overall effi- 
ciency of acceleration has been found to increase with 
the relative energy fraction carried by ions in the wind. 
The authors also show that the spectrum of non-thermal 
particles varies. 

The presence of the so-called wisps in the Crab neb- 
ula has been used to argue for the presence of ions in 
the wind (see e.g. ^Sl)- the compression caused by the 
ions leads to increased magnetic field and correspond- 
ingly more intense synchrotron radiation. Wisp-like fea- 
tures have been observed from PSR B1509-58 17, how- 
ever, the predicted time variation for the northern arc 
was not found in later observations ^3]. The Vela PWN 
shows variable features along its jet which have been in- 
terpreted to be the result of a kink instability . The 
ring like features in the Vela PWN have not been seen 
to vary with time even though the expected time scale 
should be similar to the one observed from the Crab neb- 
ula. 

Besides the existence of wisp-like structures near the 
wind shock, ions are expected to leave other observable 
signatures. The downstream energy distribution of the 
ions follows a relativistic Maxwellian distribution with 
some modifications due to the energy loss of ion energy 
transfered to the non-thermal tail of accelerated pairs |H] . 
The temperature of the distribution is close to jniiC^. 
At some distance to the shock, the ions will move diffu- 
sively outwards, loosing energy by adiabatic expansion. 
Gamma-rays and neutrinos will be produced predomi- 
nantly in inelastic scattering of the nucleons on the am- 
bient medium. Depending on the diffusion coefficient in 
the PWN, a large fraction of the particle energy can be 
converted into gamma-rays and neutrinos. 
Here, we consider specifically the Vela X PWN, which 
has recently been detected to emit VHE gamma-rays 
|17| . In the final section, we will discuss similarities and 
differences to the other "TeV Plerions" . For a more gen- 
eral review of gamma-ray production in PWN see jl8l in 
these proceedings. 



2 Vela pulsar and Vela X 

The Vela X region was initially discovered as an ex- 
tended radio emitting region , south of the Vela Pul- 
sar PSR B0833-45 located in the center of the Vela su- 
pernova remnant. The distance to the Vela pulsar is well 
known from parallax measurements to be 290 pc j20) . 
The Vela Pulsar spins with a period of 89.3 ms and slows 
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Fig. 1 In grey scale (see the online proceedings for a color 
version): ASCA 2-10 keV mosaic showing the flux in bins of 
0.25 arc min^ size. The contours represent the HESS signif- 
icance map with the contours starting at 5 cr above the 
background noise and incrementing by one standard devia- 
tion. 

down with P — 1.247 x lO^^'^ s/s which gives a current 
spin down luminosity of _E = 7 x 10'^^ erg/s and an age 
of i «11 kyrs m]. 

X-ray emission from the Vela X region was discovered 
with the ROSAT X-ray telescope |221 and initially as- 
sumed to be a "jet" from the pulsar. The spectral range 
of the ROSAT PSPC instrument was not sufficient to 
clearly identify a non-thermal tail in the observed en- 
ergy spectrum. Later high resolution measurements of 
the Vela pulsar with Chandra revealed a compact X-ray 
nebula with a double torus structure |23 . The morphol- 
ogy of the Vela compact nebula and its orientation sug- 
gests that the Vela X region is not a "jet" but rather the 
result of the interaction of the middle aged pulsar wind 
system with the reverse shock of the supernova shock 
wave. Hydrodynamic simulations of the interaction of an 
expanding PWN inside the shell remnant indicate that 
for middle aged PWN, the reverse shock will start to 
interact with the PWN leading ultimately to a compres- 
sion of the PWN ,24j|. For the specific case of asymme- 
tries in the density of the ambient medium into which 
the external shock of the supernova blastwave expands, 
the reverse shock will start to interact earlier with the 
expanding PWN along the direction where the ambient 
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Table 1 ASCA pointings used to generate the mosaic shown 
in Fig. 1 



ID Obs. date exposure 

[ksec] 



23043000 


04/15/1995 


34 


23043010 


04/15/1995 


34 


25038000 


12/02/1997 


76 


50021000 


05/12/1993 


20 


50021010 


06/26/1993 


20 


50021020 


07/14/1993 


26 


50021030 


10/08/1993 


30 



medium is more dense. The Vela X PWN shows indeed 
such asymmetry which could be the result of the inter- 
action of the evolved PWN with an asymmetric reverse 
shock. 



3 Non-thermal X-rays from Vela X: New results 
from ASCA 

The X-ray observations of the Vela X region with ROSAT 
are not conclusive on the existence of a non-thermal 
power-law component in the energy spectrumThe ASCA 
satellite with its four X-ray telescopes equipped with 
two Gas Imaging Spectrometers (GIS, and two Solid 
State Imaging Spectrographs (SIS, TP) is ideally suited 
to image an extended region like Vela X at energies 
between 2-10 keV. Given the size of the Vela X re- 
gion, we do not consider the SIS data with the smaller 
field of view. The GIS data were screened following the 
standard screening criteria. For the 4 early pointings 
(50021000, 50021010, 50021020, 50021030) no rise-time 
selection (Ohashi et al. 1996) was possible. These point- 
ings cover mainly the region at the northern end of Vela 
X. See also Table 1 for an overview of the observation 
number, observation date, and exposure (combined GIS2 
& GIS3). Since we are interested in the morphology of 
the Vela X region in non-thermal X-rays , we consider 
for the image only events with energies exceeding 2 keV. 
The data from the two GIS detectors have been added to 
increase the statistics. The particle background was esti- 
mated from Earth night sky observations and subtracted 
off the skymap and the resulting excess map with 0.5 arc 
min bins was divided by the exposure map. The analysis 
is similar to the one described in |28| . 
The resulting flux image is shown in Fig.^in grey scale. 
Note, the surface brightness is rather faint on the level of 
5 X 10~^ counts/ (s cm^ arc min^). The overlaid contours 
are from the excess map of the VHE gamma-ray source 

ini 

The X-ray morphology observed in the ASCA image is 
very similar to the ROSAT picture [221 and shows bright 
emission centered on the position of the pulsar with a 
narrow extension to the south and a re-brightening at 
the southern end of Vela X. A bright feature appears to 



the south cast of the Vela X region which is marginally 
significant and unfortunately located at the edge of the 
field of view. More observations of this region at X-ray 
energies could be of interest as there appears to be an in- 
crease of the VHE signal towards that same unexplored 
region. 

A more quantitative study of the X-ray morphology at 
different energies is beyond the scope of this paper. How- 
ever it is noteworthy, that the size of the X-ray emitting 
region at energies of 2-10 keV appears to be smaller 
than the observed size in the soft energy band as e.g. 
the ROSAT PSPC data indicate. 

In order to investigate possible spectral variations, we 
have sub-divided the extension of Vela X in three re- 
gions excluding the bright feature related to the com- 
pact nebula in the north by excising a 6 arc min radius 
region centered on the Vela Pulsar which includes the 
compact X-ray nebula [SHI- The energy spectra of all 
three regions are compatible with a mixture of a ther- 
mal component and a power-law with a photon index 
of 2 and varying fluxes. The southern tip has also been 
observed with XMM-Newton. The ASCA energy spec- 
trum of the southern tip has been cross-checked with 
the XMM-Newton spectra which show good agreement 
as both observations are well fit by a power-law compo- 
nent to be present up to 8 keV. This is consistent with 
the result of a combined ROSAT and ASCA SIS analysis 
|5n) of the southern tip region. 

In order to compare the observed X-ray emission and the 
VHE emission, an energy spectrum covering the entire 
length of the X-ray emitting region excluding a 6 arc min 
region centered on the Vela pulsar has been extracted. 
For the background estimate, the dim region to the east 
of the Vela X has been used. The resulting energy spec- 
trum is well fit by a power-law with photon index w 2 
and is shown in Fig. O together with the ROSAT flux 
and the BeppoSAX/XMM-Newton combined analysis of 
the compact nebula emission j29j. 



4 Nucleonic gamma-rays from Vela X 

4.1 Observations 

The first clear indication of a gamma-ray signal from 
the Vela X region was found with the H.E.S.S. tele- 
scopes with a luminosity in the oberved energy range 
TcV ~ 10^'^ erg/s 521 at the distance of 290 pc. 
The spatial extension of the TeV plerion is slightly larger 
than the size of the X-ray emitting region and extends 
for a full width at half maximum (FWHM) of 5.7 pc 
along the major axis and for a FWHM of 4.3 pc along 
the minor axis (again at a distance of 290 pc). The ob- 
served energy spectrum can be described by a power-law 
with a cut-off dN/dE = No{E/l TeV)"-^ • exp{-E/Ec) 
with the best fit parameters F = 1.45 ± 0.09stat ± 0.2sys 
and Ec = 13.8 ± 2.3,tat ± ^-hys TeV. 
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E[TeV] 

Fig. 2 A comparison of tlie measured VHE gamma-ray spec- 
trum |17| with tlie model fit. The dashed and dot-dashed 
curves show the expected neutrino spectra. 
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Fig. 3 Broad band spectral energy distribution from 
the Vela X region including the combined XMM- 
Newton/BeppoSAX result on the compact nebula centered 
on the Vela Pulsar ("15' X-ray nebula"), the ROSAT result 
(open triangle), the ASCA spectrum excluding the compact 
nebula ("60' TeV plerion"), and the corresponding INTE- 
GRAL upper limit from the same region (99 % c.l.) 



injection event and therefore, the particle distribution is 
assumed to follow a relativistic Maxwellian distribution 
■with a temperature /csT = jmiC^. 

When considering the relevant time-scales for the evo- 
lution of the particle distribution in comparison with 
the age of the Vela Pulsar, we assume Bohm-type dif- 
fusion and energy loss for the nucleons in inelastic scat- 
tering with the ambient medium with a density of n = 
0.6 cm^"^. This a lower limit to the overall density of this 
region derived from the emission measure of the plasma 
in X-rays. In principle, the total target density could be 
higher than this value in case of the presence of cold gas. 
The diffusion time idiff = 7300 yrsZ-iry(i?4)^B_5£;j"oJ) 
is shorter than the age of the system and therefore, es- 
cape losses are negligible. The magnetic field of 10 /uG is 
a reasonable value for the Vela X region (see also next 
section). More detailed calculation [2] indicate that the 
total energy in protons required to match the observed 
gamma-ray flux is Wp « lO''^ erg. 

A fit of the energy spectrum of gamma-rays from proton- 
proton interaction assuming a relativistic Maxwellian dis- 
tribution is shown in Fig. El The 7r° decay spectrum 
has been calculated using the parametrization given in 
122] . The observed energy spectrum is well described for 
E = 80 TeV as the energy of the Maxwellian which trans- 
lates into 7 « 8 X 10'*. 

A characteristic signature of the nuclconic production 
mechanism is the production of neutrinos. Using the cor- 
responding parametrization for the production of muon 
and electron flavor neutrinos [^2], we calculate the dif- 
ferential energy spectra for the two flavors also shown in 

Fig.m 

Given the high energy cut-off and the reasonable angular 
extension of the Vela X region, we consider Vela X among 
the best candidates for a detection as a neutrino source 
with the future Neutrino telescope in the mediterranean 
sea. 



4.3 X-ray emission from secondary electrons 



4.2 Model 

The observed gamma-ray spectrum can be modeled un- 
der the assumption that the pulsar wind carries a large 
fraction of the spin down energy in the form of nucle- 
ons that expand diffusively to fill the PWN. The pulsar 
releases most of its energy during the first few hundred 
years, when the spin down proceeds most rapidly. While 
the pulsar spins down, ^opon decreases and therefore it is 
natural to assume that the corresponding Lorentz factor 
of the particle wind decreases as well, even though we 
do not have as yet direct evidence for this evolutionary 
effect. 

Given the short time in which the particles are injected, 
it is reasonable to assume that we have an instantaneous 



Due to the interaction of nuclei with ambient medium, 
neutral and charged mesons are produced. While the 
neutral mesons (primarily tt") decay to produce 7-rays, 
charged pions decay into electrons and positrons. These 
secondary electrons then in turn radiate mainly via syn- 
chrotron, inverse Compton, and Bremsstrahlung. The 
energy loss time of the electrons depends on the energy 
density of the background radiation field (magnetic field 
strength and soft photon density) and the particle num- 
ber density. In the case of Vela X, synchrotron cooling 
dominates very likely with a cooling time expressed as a 

function of the characteristic energy of synchrotron radi- 

3/2 1/2 

ation emitted: ti/2 = 1.2 kyrsi?_g Ef.^y . The magnetic 
field strength in the Vela X region is difficult to estimate. 
The equipartition field derived from radio measurements 
of narrow structures can be estimated to be 20 — 50 /iG. 
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However, when taking the X-ray flux in the considered 
region, the equipartition field is B^q = 4 /iG. For such 
a small magnetic field however, electrons are not effi- 
ciently confined in the Vela X region and should have 
escaped already and would fill a larger volume which in 
turn should produce a considerably larger TeV plerion. 
Therefore, a magnetic field around 10 /zG is considered 
to be a realistic parameter to match the observed size of 
the TeV plerion. 

The synchrotron emission from secondary electrons is 
calculated taking into account radiative cooling but ne- 
glecting escape losses (see above). The resulting syn- 
chrotron spectrum assuming a 10 /iG magnetic field after 
11000 yrs is shown in Fig|3| together with the relevant 
measurements. Clearly, the secondary electrons' contri- 
bution to the overall spectral energy distribution is neg- 
ligible. In principle, the synchrotron component would 
dominate in the optical. However, the sensitivity for ex- 
tended emission from the Vela X region in the optical is 
at least one order of magnitude above the expectation 

m 



4.4 X-ray emission from primary electrons 

In the framework of ion resonant acceleration, pairs are 
accelerated to follow a power-law which extends up to a 
maximum energy given by the energy of the upstream 
ions. With the observed gamma-ray spectrum, the up- 
stream ion energy is constrained to be of the order of 
100 TeV. Given the energy losses of the downstream ions, 
the upstream value (jmiC^) is higher than the temper- 
ature of the Maxwellian downstream distribution which 
is subsequently modified by adiabatic losses. We assume 
for simplicity that the pulsar is currently injecting a 
wind with a Lorentz factor of « 10^ infered from the 
fit of the observed gamma-ray spectrum neglecting time- 
dependent effects. 

The primary electrons accelerated at the shock will there- 
fore reach a maximum energy of E'max ~ 100 A/Z TeV 
with Z the charge and A the mass number of the nuclei in 
the wind |Hj . Given the magnetic field in the downstream 
region (see below) , this value limits the maximum energy 
of emitted synchrotron radiation. 

The magnetic field at the shock can be calculated in the 
standard MHD picture of PWN (e.g. |33]) by consider- 
ing the distance of the shock to the pulsar, the total 
extent, and luminosity of the nebula. The most recent 
estimates range from B{6s = 33") = 72 /iG j25l as- 
suming the position of the shock at an angular separa- 
tion of 33 arc sec while using the value of 21 arc sec 
obtained by an elaborate fitting method of the torus 
IHE] increases B{9s = 21") = 113 fiG. This value is 
obtained assuming a magnetization of u = 0.1 which 
appears consistent with MHD simulations modeling the 
size of the Vela nebula |37,. Considering the shock com- 
pression the downstream magnetic field can reach values 



of B„,ax « 2 5(0,) « 200 fiG. 

Correspondingly, the emitted synchrotron radiation of 
electrons is given by fmax = 50 keV{B-4){A/Z)'^ E^qq. 
The detection of unpulsed X-ray emission up to 200 keV 
with INTEGRAL from the Vela pulsar/PWN 38 indi- 
cates that possibly the nucleons would have to be only 
partially ionised to increase (A/Z) and therefore the max- 
imum energy of the pairs. Partially ionised nuclei have 
already been suggested in [HI in order to avoid exceed- 
ing the Goldreich-Julien current. It is quite interesting 
to note that most of the pulsars powering the TeV pleri- 
ons have been also detected as INTEGRAL hard X-ray 
sources (see Hoffmann et al. these proceedings). More 
data on hard X-ray emission from TeV plerions is of great 
importance to understand the in situ acceleration at the 
pulsar wind shock and its maximum energy. 
The Vela X region has not been detected to emitt X- 
rays beyond « 8 keV. Taking all available data with 
INTEGRAL results in a conservative upper limit on the 
energy flux in the 20-60 keV band as shown in Fig. |21 
Comparing the upper limit with the spectrum obtained 
with ASCA it is evident, that a cut-off in the X-ray spec- 
trum with an energy « 10 keV is required not to violate 
the upper limit. 

The origin of the X-rays is very likely synchrotron emis- 
sion from the primary electrons that are accelerated at 
the shock and show a radiative cooling break which moves 
from a few keV to below keV energies when extracting 
the energy spectra at increasing angular separation from 
the Vela pulsar (3^1. This is consistent when consider- 
ing the Vela X region, where it turns out to be below 
keV energies. Under the assumption of a constant injec- 
tion rate with a power-law with index 2, the synchrotron 
spectrum for a 10 /xG field after 11000 yrs is shown in 
FigEl The total injected energy in electrons amounts 
to We ~ 10"^ B^^ erg between E^in = 0.01 TeV and 
^'max — 200 TeV with the minimum energy chosen close 
to jmeC^ = 0.5 TeV. 



5 Other TeV plerions 

Besides the Vela X TeV plerion, at least two more mid- 
dle aged PWN have been observed to be TeV pleri- 
ons: G18.0-0.7 3& and G313.3-f0.6 40; (the northern 
wing of the "Kookaburra" ) . It is interesting to point out 
similarities and distinct differences of these objects: All 
three objects are spatially resolved to be extended at 
VHE gamma-rays. While the Vela X region measures 
« 5 pc across, G18.0-0.7 (25.3 pc at 4 kpc distance) 
and G313.3-I-0.6 (12.1 pc at 5.6 kpc distance) fill a sub- 
stantially larger volume. The differences in size can be 
at least partially attributed to the age: the Vela pulsar 
{t = 11 kyrs) is nominally the youngest of the three 
objects (G18. 0-0.7 is 21.5 kyrs and G313.3+0.6 13 kyrs 
old). However, there are indications from the very small 
braking index of the Vela pulsar that the actual age of 
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Vela is closer to 50 kyrs which would be more consistent 
with the rather cold neutron star surface. 
When comparing the spin down power of the pulsar with 
the radiative power (either in X-rays or gamma-rays), 
it is interesting to note that while in Vela, the radia- 
tive power accounts for less than 1 per cent of the spin 
down power, it accounts for a few per cent in the other 
objects. Finally, we consider the ratio of the published 
extension (gg for G18.0-0.7 and for G313. 3+0.6) 
aof the X-ray nebula rx to the radius of the TeV ple- 
rion TTeY- for Vela X, this ratio is fairly close to 1 while 
it is only rx/rxeV « 0.1 for G18.0-0.7 and G313.3-h0.6. 
It is beyond the scope of this paper to investigate the 
deeper reason for this large difference, but it is clear 
that the spatially resolved gamma-ray spectroscopy with 
H.E.S.S. combined with deep X-ray observations cover- 
ing the same angular region will greatly improve the un- 
derstanding of the underlying physical parameters. 
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